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I.  INTRODUCTION 


A  radio  frequency  (if)  discharge  dissodator  is  die  most  frequently  used  source  of  atomic  hydro¬ 
gen.  Hydrogen  masers,  the  time  and  frequency  standards  of  choice  whenever  the  highest  possible 
frequency  stability  is  required,  use  an  if  dissodator  followed  by  a  hexapole  or  quadrupole  mag¬ 
net  to  focus  a  beam  of  state-sdected  hydrogen  atoms  at  die  entrance  orifice  of  die  maser  bulb. 
Similar  combinations  of  if  discharge  dissodators  and  state-selecting,  focusing  magnets  have  been 
used  to  produce  beams  of  spin-polarized  hydrogen  atoms  used  in  atomic,  nuclear,  and  high- 
energy  physics  experiments.2’3  Since  atomic  trajectories  within  a  hexapole  or  quadrupole  mag¬ 
net  will  depend  strongly  on  atomic  velodties,  such  magnets  also  act  as  fairly  narrow-band  velodty 
filters.4  Efficient  use  of  die  state-sdected  atomic  hydrogen  source,  which  becomes  a  paramount 
consideration  when  designing  hydrogen  masers  qualified  for  space  flight,  requires  a  good  match 
between  the  vdodty  distribution  of  die  hydrogen  atoms  exiting  die  dissodator  and  die  velocity 
bandpass  of  die  state-seteedng  magnet 

In  many  instances,  die  tadt  assumption  is  made  that  die  atoms  exiting  die  dissodator  will  be  in 
thermal  equilibrium  with  the  dissodator  wall;  this  assumption  is  not  necessarily  valid. 

Dissodadon  of  molecular  hydrogen  in  an  rf  discharge  taku  {dace  after  spin-exchange  collisions 
with  energetic  discharge  electrons.3  The  ground  state  of  the  H2  molecule  has  a  configura¬ 
tion;  die  lowest  energy  triplet  stale  is  die  repulsive  3£Q,  which  has  the  same  asymptotic  dissoda¬ 
don  limit  as  die  ground  state  at  about  4  J  eV,  and  which  also  will  be  the  endpoint  of  any  radiative 
cascade  starting  after  exdtation  of  a  higher  triplet  level.  Thus,  any  spin-exchange  collision  will 
end  up  dissodating  the  molecule.  Since  die  threshold  energy  for  exdtation  of  the  3I„  state  is  8.8 
eV,6  after  dissociation,  each  one  of  die  hydrogen  atoms  will  carry  away  more  than  2  eV  of  excess 
energy.  Depending  on  dissodator  geometry  and  gas  density,  the  hydrogen  atoms  may  or  may 
not  undergo  enough  gas  and  wall  collisions  to  thermalize  frilly  before  they  exit 

Relatively  little  is  actually  known  about  die  velodty  distribution  of  hydrogen  atoms  effusing  out 
of  rf  discharge  dissodators.  Walraven  and  Silvers '  have  studied  the  characteristics  of  a  beam  of 
hydrogen  atoms  produced  in  a  microwave  discharge  followed  by  a  cryogenic  thermal  accom¬ 
modator  and  found  die  velodty  distribution  to  be  Maxwellian  at  the  accommodator  temperature. 
Hershcovitch  et  aL 8  studied  a  similar  dissodator- accommodator  combination  operated  at  higher 
pressures  and  flow  rates  and  measured  vdodty  distributions  in  reasonable  agreement  with  the 
results  of  supersonic  flow-regime  calculations.  These  results  show  that  accommodators  can  ther¬ 
malize  the  hydrogen  beam,  but  they  are  not  directly  applicable  to  die  problem  under  study  since 
thermal  accommodators  are  not  commonly  used  in  hydrogen  masers  or  polarized  hydrogen 
beam  sources.  Miner9  investigated  a  helium-cooled  microwave  discharge  dissodator  at  relatively 
Ugh  pressure  (3.3  tore)  and  measured  a  velodty  distribution  that  was  slightly  narrower  than  a 
beam-Maxwellian  at  the  estimated  discharge  temperature. 

We  have  determined  velodty  distributions  of  atomic-hydrogen  beams  effusing  out  of  an  rf  dis¬ 
charge  dissodator  having  a  geometry  and  operating  parameters  similar  to  those  of  a  maser  disso¬ 
dator  by  analyzing  the  deflection  of  the  atomic  trajectories  in  an  inhomogeneous  magnetic  field. 
These  distributions  are  presented  and  compared  with  the  results  of  Monte  Carlo  simulations  of  the 
atomic  hydrogen  thermalization  process  within  the  dissodator  bulb. 
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II.  EXPERIMENTAL  TECHNIQUE 


Atomic  hydrogen  velocity  distributions  have  been  determined  using  a  magnetic  deflection  tech¬ 
nique.  Figure  1  shows  schematically  our  experimental  arrangement.  Hydrogen  is  fed  through  a 
standard  temperature-controlled  Pd-Ag  leak*0  into  a  cylindrical  double-walled  Pyrex  bulb  15 -cm 
long  with  a  1.9 -cm  internal  diameter  ami  a  3.8-cm  external  diameter.  Compressed  air  flows 
between  the  internal  and  external  walls  to  provide  cooling,  and  rf  power  is  coupled  inductively  to 
the  discharge  by  an  external  25-turn  coil.  The  hydrogen  beam  exits  the  dissodator  through  a 
0.025 -cm-wide,  0.1-cm-high  slit,  is  collimated  by  a  second  0.025-cm-wide  slit  set  at  d  =  63.7  cm 
away  from  die  dissodator  slit,  and  then  travels  between  die  polepieces  of  an  L  =  11.4-cm-long, 
inhomogeneous-field  dectromagnet,  configured  in  the  "two-wire"  geometry  described  by  Rabi  et 
al11.  After  additional  travel  through  a  D  *  71.3  cm  drift  tube,  die  beam  is  detected  by  an  rf 
quadrupole  mass  analyzer.  The  dissodator  is  attached  to  the  rest  of  the  apparatus  by  a  flexible 
vacuum  bellows,  and  can  be  displaced  transversely  by  micrometer  screws,  allowing  measurement 
of  the  angular  distribution  of  atoms  deflected  by  die  inhomogeneous  magnetic  field. 

The  atomic  traction  in  the  hydrogen  beam,  <p,  was  measured  using  the  quadrupole  mass  analyzer. 
The  results  were  corrected  for  the  H+  fraction  when  detecting  H2;  typically,  values  between  0.65 
and  almost  1.00  were  obtained,  depending  on  total  dissodator  pressure  and  rf  power  coupled  to 
the  discharge.  The  dissodator  pressure  was  measured  using  a  capadtance  manometer,  and 
injected  and  reflected  rf  power  were  monitored  continuously  as  indicators  of  die  steady-state 
character  of  die  discharge. 


-* — d — L — - D - H 


DISSOCIATOR  MAGNET  MASS  FILTER 

AND  DETECTOR 

Figure  1.  Schematic  representation  of  the  experimental  arrangement 
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Ail  atom  having  magnetic  moment,  p,  immersed  in  a  magnetic  field,  B,  has  energy,  E  =  -p«B.  The 
force  acting  on  the  atom  is  F  =  grad(p*B)  *  iiefigrad(B),  where  peff  is  foe  component  of  the 
atomic  magnetic  moment  along  foe  direction  of  foe  field.  In  general,  peff  will  depend  on  foe 
atomic  ground-state  quantum  numbers  and  the  field  strength.  For  an  atom  having  total  electron 
angular  momentum  and  nuclear  spin  quantum  numbers  J  *  1/2  and  I,  respectively,  and  total 
azimuthal  quantum  number,  M,  foe  effective  magnetic  moment,  is  given  by  Breit-Rabi's 
formula12. 


H«ff=± 


e+2M/(2I+l) 

[1 + 4Me  /  (21 + 1) + e2  ]y2  ^°’ 


O) 


where  foe  (+)  sign  corresponds  to  foe  total  angular  momentum  quantum  number,  F  =  1-1/2,  and 
foe  (-)  sign  to  F=  1+1/2.  The  parameter  e  =  gp^B/AW,  where  g  is  foe  atomic  Lande  factor,  AW  is 
foe  atomic  hyperfine  splitting,  and  pg  is  toe  Bohr  magneton. 

Figure  2  shows  foe  cross  section  of  foe  "two-wire"  magnet;  2a  is  foe  separation  between  foe 
"equivalent  wires."13  If  a  narrow  atomic  beam  enters  foe  magnet  at  x  =  1.2a  (as  determined  by  a 
collimating  slit),  and  the  atomic  transversal  displacements  within  foe  magnet  are  small,  foe  force 
acting  on  foe  atoms  remains  approximately  constant,14  and  so  does  foe  acceleration: 


Figure  2.  Cross  section  of  foe  "two-wire*  magnet  polepieces. 
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,.09*1^ 
a  m 

where  B  is  the  magnetic  field  strength  at  the  atomic  beam,  m  is  the  mass  of  the  atom,  and  i  a  unit 
vector  along  die  x  (deflection)  axis.  Under  these  conditions,  the  trajectory  equations  for  the 
atoms  traveling  from  source  to  detector  can  be  solved  easily.  With  no  magnetic  field  applied,  an 
atom  of  effective  magnetic  moment  and  speed  v  leaving  the  source  slit  at  x§  will  be  detected 
at  xD  *  -xs(L+D)/d.  If  now  a  magnetic  field  B  is  applied,  and  the  detector  is  not  moved,  an 
identical  atom  will  have  to  leave  the  plane  of  the  source  slit  at  x§',  given  by 

*Xs+^t(,+l?d)  (3) 

in  order  to  be  detected  at  xD.  If  fQ(x§)  is  the  distribution  of  detected  atoms  as  a  function  of 
source-slit  position  at  zero  field,  and  g(v)  is  the  speed  distribution  of  atoms  leaving  the  source  slit, 
the  distribution  of  detected  atoms  at  non-zero  field  vs.  source-slit  position,  f(x$),  will  be  given  by 


f(xs)=  XJT  fo(*S  “kfm  /  v^gCvjdv,  (4) 

F.M 


where 

Both  fQ(xg)  and  f(xs)  can  be  measured,  and  then  Eq.  (4)  can  be  used  to  verify  whether  a  model 
velocity  distribution,  g(v),  is  consistent  with  the  measured  deflected  atom  distribution. 
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III.  CALIBRATION  AND  VERIFICATION  EXPERIMENT 


A  preliminary  experiment  involving  the  use  of  a  rubidium  atomic  beam  was  performed  in  order 
to  calibrate  the  electromagnet  and  test  the  use  of  Ea.  (4)  to  verify  a  model  velocity  distribution. 
A  two-chambered  effusive  oven  of  standard  design‘d  was  installed  at  the  end  of  the  source  bel¬ 
lows,  where  the  rf  dissociate  would  normally  be  located. 

As  shown  by  Eq.  (1),  an  atom  will  have  =  0  for  B  such  that  £  =  -2M/(2I+1).  At  field  intensi¬ 
ties  determined  by 


-M  AW 

21+1  Ho 


(6) 


atoms  with  M  <  0  will  not  be  deflected  by  die  inhomogeneous  field  in  the  magnet  Natural  rubid¬ 
ium  is  a  mixture  of  72.2%  ®%b  (I  =  5/2)  and  27.8%  ^Rb  (I  =  3/2).  *^Rb  has  zero  effective 
moments  at  B  =  361G  (M  =  -1)  and  B  =  722G  (M  =  -2),  while  has  zero  effective  moments 
at  B  =  1221G  (M  =  -1).  Figure  3  shows  the  detected  atomic-beam  signal  as  a  function  of  mag¬ 
net-coil  current,  with  the  rubidium  beam  source  on  die  apparatus  axis.  At  die  field  values  listed 


Figure  3.  Rubidium  atomic-beam  signal  (oven  and  detector  oo  axis)  vs. 

magnet-coil  current  Arrow  shows  features  due  to  zero  effective 
moment  states.  Dasbed  lines  were  obtained  by  subtracting  the 
local  constant-slope  background  from  each  feature.  A:  Rbg,,  M 
*-2;  B:  Rbg7,  M  -1. 
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above,  atoms  in  the  Pgg-  =  0  states  will  travel  through  the  "two-wire"  magnet  without  deflection, 
leading  to  an  increased  detector  signal.  The  effect  of  the  zero  effective  moment  states  can  be 
made  more  apparent  in  Figure  3  by  subtracting  from  each  one  of  the  marked  features  of  the  local 
sloping  background.  The  peaks  labeled  A  and  B  can  be  identified  with  the  zero  effective 
moments  at  722G  and  1221G,  respectively,  yielding  a  calibration  of  66.3G/A  for  the  "two-wire" 
magnet. 

The  rubidium  atomic  beam-speed  distribution  was  assumed  to  be  beam-Maxwellian, 


g(v)=  (2/p)  (v/P)3  exp[-(v/p)2],  (7) 

where  P  =  ^kT/rn)1/2,  and  T  =  484  K,  the  absolute  oven  temperature.  We  measured  f0(xs)  and 
f(xs),  and  then  used  Eq.  (4),  with  the  proper  statistical  weight  for  each  rubidium  isotope,  to  calcu¬ 
late  the  expected  f(x§); 


f(x<j)  —  0.722fgj(Xg)  +  0.278fg7(xs). 


(8) 


Figure  4  shows  measured  and  calculated  deflected-atom  distributions  at  B  =  106 1G.  For  xs  <  0 
the  agreement  between  measured  and  calculated  f(x§)  is  excellent;  these  results,  showing  that  the 
rubidium  beam-speed  distribution  was  indeed  a  source-temperature  beam-Maxwellian,  validate 
our  experimental  technique.  The  agreement  is  poorer  for  x§  >  0;  this  is  the  convex  polepiece 
side  of  the  magnet,  and  atoms  may  be  traveling  too  close  to  it  Thus,  only  the  xs  <  0  side  of  the 
deflection  pattern  was  considered  when  evaluating  model-speed  distributions  for  the  atomic- 


Figure  4.  Distribution  of  deflected  rubidium  atoms  at  B  =  1061  Gauss.  Full 
line:  measured  distribution.  Dashed  line:  distribution  calculated 
assuming  a  beam-Maxwellian  atomic-speed  distribution. 
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IV.  ATOMIC  HYDROGEN  DEFLECTION  MEASUREMENTS 


Since  the  focus  of  this  study  was  on  the  degree  of  thermalization  of  the  atomic-hydrogen  beam 
effusing  out  of  die  dissodator,  we  used  our  magnetic  deflection  technique  to  investigate  die 
dependence  of  the  atomic-hydrogen  speed  distributions  on  dissociate  pressure.  Deflection 
profiles  f(x§)  were  measured  at  B  =  995 G  (corresponding  to  15  A  excitation  current)  for  total 
dissociate  pressures  between  28  and  340  mtorr.  Undeflected  (zero-field)  beam  profiles,  fg(xs), 
and  atomic-hydrogen  beam  fractions,  9,  were  also  measured  in  each  case. 

Unde  free  molecular  flow  conditions,  fQ(x§)  should  be  fully  determined  by  the  source  and  col¬ 
limator  slits,  and  be  independent  of  source  pressure.  This  was  indeed  the  case  for  the  range  of 
dissociate  pressures  we  have  explored,  as  shown  in  Figure  5.  The  line  in  the  figure  shows  %(x§) 
calculated  for  our  geometry  assuming  a  trapezoidal  beam  intensity  distribution*6  folded  with  a 
detector  of  finite  width;  the  calculation  contains  no  free  parameters  and  is  in  excellent  agreement 
with  the  measurements.  The  slight  discrepancies  at  the  top  (scans  proceeding  from  left  to  right) 


Figure  5.  Atomic-hydrogen  beam  profile  (zero  deflecting  magnetic  field). 

Squares:  measurements  at  3 19  mtorr.  Triangles:  measurements  at 
215  mtorr.  Hexagons:  measurements  at  58  mtorr.  Full  line:  cal¬ 
culated  beam  profile. 


are  due  to  detection  system  response/recovery  times  which  were  comparable  to  source  point- 
dwell-times.  We  believe  the  calculated  fQ(x§)  to  be  a  better  description  of  the  undeflected  beam 
than  any  of  our  measurements;  thus,  it  was  used  with  Eq.  (4)  to  calculate  deflected  beam 
distributions. 

Deflected  beam  profiles  were  calculated  using 

f(xs)  -  <pft(xs)  +  (l-q»f0(xs)  ,  (9) 

where  the  first  term  describes  the  deflected  atomic-hydrogen  beam  profile,  with  fj  (x$)  computed 
using  Eq.  (4),  and  the  second  term  represents  the  contribution  of  the  H2  fraction  of  the  beam, 
which  is  not  deflected  by  the  magnet  We  performed  a  non-linear,  least-squares  fit  to  the  x§  <  0 
deflection  data  using  Eq.  (9),  with  a  beam-Maxwellian,  g(v),  Eq.  (7),  and  0  as  a  free  parameter. 
Figure  6(a)  shows  the  results  for  p  s  319  ratorr.  The  left  lobe  of  the  best  fitting  f(xs),  obtained 
for  0  *  2.29  x  103  m/s  and  indicated  by  the  full  line,  is  in  quite  good  agreement  with  the  mea¬ 
sured  deflected  beam  data,  indicating  that  therm  alizahon  of  the  hydrogen  atoms  within  the  disso- 
ciator  bulb  is  essentially  completed.  The  equilibrium  temperature  derived  from  the  best  value  of 
0  is  318  K.  The  results  of  the  fits  at  lower  dissotiator  pressures,  on  the  other  hand,  strongly  sug¬ 
gest  incomplete  foermalization.  Figure  6(b)  shows  measured  and  best-fitting  deflected-beam 
profiles  for  p  s  28  mtorr.  Inspection  of  the  left  lobe  of  the  beam  profile  shows  that  the  f(xs)  cal¬ 
culated  using  the  beam-Maxwellian,  g(v),  indicated  by  file  full  line,  predicts  too  many  atoms 
undergoing  very  small  as  well  as  very  large  deflections.  Since  atomic  deflections  Ax§  «  l/v2,  this 
suggests  that  the  low-  and  high-speed  tails  of  the  beam-Maxwellian  speed  distribution  are  sub¬ 
stantially  higher  than  those  of  the  actual  atomic-hydrogen  beam-speed  distribution  at  low  disso- 
ciator  pressure.  Similar  results  are  shown  in  Figure  6(c)  for  p  =  86  mtorr. 

In  order  to  evaluate  the  differences  between  the  actual  atomic-hydrogen  beam-speed  distribution 
and  the  beam-Maxwellian  distribution,  which  would  be  characteristic  of  a  frilly  thermalized 
atomic-hydrogen  beam,  we  used  the  model  distribution 


g(v)  =  CN  v3  exp[-(v-vQ)2/2cr2], 


(10) 


where  Vq  and  <r  are  treated  as  free  parameters  in  a  non-linear,  least-squares  fitting  procedure. 
is  a  normalization  factor,  and  g(v)  becomes  beam-Maxwellian  in  the  limit  vq  =  0,  in  which  case  0 
=  v2o.  Best-fitting  deflected-beam  profiles,  f(xs),  obtained  using  that  model  distribution  are 
shown  in  Figure  6  by  the  dashed  lines.  The  results  at  319  mtorr  are  comparable  to  those  obtained 
with  the  beam-Maxwellian  distribution,  but  at  the  lower  pressures,  28  and  86  mtorr,  they  are  in 
much  better  agreement  with  the  measured  deflection  profiles. 

The  model  speed  distributions  obtained  using  Eq.  (10)  with  the  best-fitting  values  of  vq  and  a  fen' 
p  a  319,  28,  and  86  mtorr  are  shown  in  Figures  7(a),  (b),  and  (c),  respectively,  together  with  the 
corresponding  best-fitting  beam-Maxwellians.  As  expected  from  the  previous  discussion,  at  the 
lower  pressures,  the  best-fitting  model  distributions  are  significantly  narrower  than  foe  beam 
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NORMALIZED  SIGNAL 


Maxwellian;  at  higher  pressures,  the  two  become  indistinguishable,  indicating  full  thermal ization 
of  the  atomic-hydrogen  beam. 


Figure  6.  Distribution  of  deflected  hydrogen 
atoms  at  B  =  995  Gauss.  Squares; 
measurements.  Full  lines: 
calculations  using  die  best-fitting 
beam-Maxwellians.  Dashed  lines: 
calculations  using  two-parameter 
model  speed  distributions,  (a)  p  = 
319  mtorr;  (b)  p  =  28  mtorr;  (c)  p 
=  86  mtorr. 


Figure  7.  Atomic-hydrogen  speed 

distributions  giving  best  fits  to  die 
magnetic-deflection  data.  Full 
lines:  beam-Maxwellian 
distributions.  Dashed  lines:  two- 
parameter  model  distributions,  (a) 
p  =  319  mtorr,  (b)  p  =  28  mtorr; 
(c)  p  =  86  mtorr. 
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V.  THERMAUZAT10N  PROCESS  SIMULATIONS 


The  low-pressure  results  presented  in  the  preceding  section  confront  us  with  an  apparent  paradox: 
since  the  atomic  beam-speed  distribution  is  non-Maxwellian,  the  hydrogen  atoms  in  the  beam 
must  leave  the  dissodator  bulb  without  undergoing  enough  collisions  to  thermalize  fully,  but,  at 
the  same  time,  they  must  undergo  enough  collisions  to  drop  from  an  initial  energy  of  more  than 
2  eV  (as  discussed  in  Section  I)  to  die  0.05  eV  characteristic  of  die  measured  speed  distributions. 
These  seemingly  contradictory  requirements  can  be  satisfied  simultaneously  in  die  presence  of  a 
process  having  a  significant  collision  cross  section  and  a  large  energy  loss  per  event;  impact 
vibrational  excitation  of  the  H2  molecules  present  in  the  dissodator  (a  loss  of  about  0.S2 
eV/collision17)  is  a  likely  candidate,  allowing  the  kinetic  energy  of  the  hydrogen  atoms  to  drop  to 
very  low  values  in  about  four  such  collisions,  well  before  a  fully  thermal  distribution  can  be 
developed.  In  order  to  provide  qualitative  support  for  this  interpretation  of  the  speed  distribu¬ 
tions  inferred  for  low  dissodator  pressures,  we  developed  a  computer  code  to  perform  Monte 
Carlo,  two-dimensional  simulations  of  the  thermalization  of  hydrogen  atoms  produced  by  disso¬ 
ciation  of  3IU  H2  molecules. 

In  our  simulations,  hydrogen  atoms  with  an  initial  kinetic  energy  of  2.15  eV  are  bom  at  random 
positions  within  a  two-dimensional  rectangular  enclosure  having  an  exit  aperture.  The  enclosure 
and  aperture  dimensions  are  those  of  our  dissodator.  The  atoms  are  followed  as  they  travel 
within  the  enclosure,  colliding  with  other  hydrogen  atoms,  hydrogen  molecules,  or  the  enclosure 
walls.  Energy  is  exchanged  in  elastic  collisions  with  hydrogen  atoms  or  molecules,  inelastic 
(rotational  or  vibrational)  collisions  with  the  molecules,  or  thermal  accommodation  at  the  enclo¬ 
sure  walls.  Each  atom  is  followed  until  it  exits  through  the  aperture  or  recombines  on  the  enclo¬ 
sure  boundary.  If  the  atom  exits  the  enclosure  within  the  atomic  beam  acceptance  angle,  the 
count  in  the  appropriate  speed  bin  in  the  atomic  beam-speed  distribution  is  incremented  by  one. 
The  simulation  stops  after  a  preset  number  of  beam  atoms  have  been  counted.  The  input 
parameters  are  the  wall  temperature,  the  total  hydrogen  pressure,  the  hydrogen  dissodation  frac¬ 
tion,  die  cross  sections  fin  (H,H)  and  (HJ^)  elastic  collisions  as  well  as  (H,H2)  vibrational  and 
rotational  exdtations,  and  the  average  thermal  accommodation  coeffident  in  collisions  with  the 
wall.  The  output  of  the  simulation  is  a  statistical  sampling  of  the  atomic-hydrogen  beam-speed 
distribution,  integrated  over  400-m/s-width  speed  bins. 

The  simulations  were  performed  fen  a  two-dimensional  gas  in  order  to  shorten  the  computation 
time  as  much  as  possible,  but  they  incorporate  all  of  the  relevant  physical  processes.  The  results 
obtained  fin  300, 30,  and  90  mtorr  are  shown  in  Figures  8(a),  (b),  and  (c),  respectively.  The 
beam-Maxwellian  speed-distribution  characteristic  of  a  two-dimensional  gas, 

g(v)=^(v/P)2expKv/0)2],  (11) 
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Figure  8.  Atomic-hydrogen  speed  distributions  for  a  two-dimensional  gas, 

integrated  over  the  width  of  a  400-m/s  speed  bin.  Squares:  Monte 
Carlo  simulation  results.  Full  lines:  beam-Maxwellians  at  tbe 
320  K  wall  temperature,  (a)  p  =  300  mtorr,  (b)  p  =  30  mtorr,  (c)  p 
=  90  JUorr. 
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calculated  at  the  dissodator  wall  temperature,  T  «  320  K,  and  integrated  over  the  400-m/s  speed- 
bin  width,  is  also  shown  in  each  figure  for  comparison  purposes.  At  the  two  lower  pressures,  the 
simulation  results  are  significantly  narrower  than  the  two-dimensional  beam  Maxwellian.  At  300 
mtorr,  file  simulation  results  are  in  excellent  agreement  with  the  two-dimensional  beam 
Maxwellian  (allowing  for  the  statistical  spread  characteristic  of  the  results  of  a  Monte  Carlo  calcu¬ 
lation),  indicating  a  completely  thermalized  atomic-hydrogen  beam.  These  results  are  fully  con¬ 
sistent  with  our  experimental  determinations  of  atomic-hydrogen  beam-speed  distributions.  If 
the  vibrational  excitation  energy  loss  mechanism  is  "turned  off"  in  the  simulation  by  assigning  a 
very  small  value  to  the  corresponding  cross  section,  very  different  results  are  obtained  for  low 
pressures,  as  illustrated  in  Figure  9.  The  flat,  fast-atom  tail,  caused  by  foe  lack  of  an  efficient 
energy  loss  mechanism,  would  show  up  in  the  magnetic-deflection  measureme  tained  for 

low  dissodator  pressures  as  a  significant  signal  excess  near  beam  axis;  no  such  s  is  observed 

in  Figures  6(b)  and  (c). 


SPEED  (M/S) 

Figure  9.  Atomic-hydrogen  speed  distributions  fcr  a  two-dimensional  gas, 

integrated  over  the  width  of  a 400-m/s  speed  bin.  Squares:  Monte 
Carlo  simulation  results  at  30  mtorr  dissodator  pressure,  with 
(HJI2)  impact  vibrational  excitation  disallowed.  Bill  line:  beam- 
MaxweQian  at  320  K  wall  temperature. 
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VI.  CONCLUSIONS 


The  results  of  this  study  show  that  the  Maxwellian  character  of  the  velocity  distribution  in  an 
atomic-hydrogen  beam  effusing  out  of  an  rf  discharge  dissociator  cannot  be  taken  for  granted. 
Full  tbermalization  will  always  be  achieved  for  high  enough  dissociator  pressures,  but  when 
operating  at  low  dissociator  pressures,  non-thermal  speed  distributions  may  obtain.  Of  course, 
what  pressure  constitutes  a  "low"  pressure  will  be  determined  in  each  case  by  the  dissociator 
geometry.  The  unavoidable  presence  of  some  molecular  hydrogen  within  the  dissociator  bulb 
makes  the  escape  of  high-energy  atoms  very  unlikely  since  impact  vibrational  excitation  of  H2 
provides  a  very  efficient  energy-loss  mechanism;  thus,  at  low  pressures,  die  atomic-speed  distri¬ 
bution  lacks  any  significant  "hot-atom"  tails  that  might  be  expected  otherwise.  The  speed  distri¬ 
butions  that  we  have  inferred  from  our  magnetic -deflection  measurements  have  averages  that  are 
quite  dose  to  thermal,  but  at  low  dissociator  pressures,  they  are  significantly  narrower  than  beam- 
Maxwellians. 
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•E13C2MOLOGY  OPERATIONS 


The  Aerospace  Corporation  ftecte  cs  aa  “architect-engineer"  for  uadcaal  security 
programs,  specialising  ia  advanced  military  space  systems.  The  Corporation's  Technology 
Operations  supports  the  effective  and  timely  development  and  operation  of  national  security 
systems  through  scientific  rcsserch  and  the  application  of  advanced  technology.  Vital  to  the 
success  of  the  Corporation  is  the  technical  staffs  wide-ranging  expertise  and  its  ability  to  stay 
abreast  of  new  technological  developments  end  progress  support  issues  associated  with  rapidly 
evolving  space  systems.  Contributing  capabilities  are  provided  by  these  individual  Technology 
Centers: 

Electrodes  Testeology  Center:  Microelectronics,  solid-state  device  physics, 

VLSI  reliability,  compound  semiconductors,  radiation  hardening,  data  storage 
techao’ogics,  infrared  dstector  devices  and  testing;  electro-optics,  quantum  electronics, 
solid-stats  lasers,  optical  propagation  and  communications;  cw  and  pulsed  chemical 
laser  development,  optical  resonators,  beam  control,  atmospheric  propagation,  and 
loser  effects  and  countermeasures;  atomic  frequency  standards,  applied  laser 
spectroscopy,  laser  chemistry,  laser  optoelectronics,  phase  conjugation  and  coherent 
imaging,  sols?  cell  physics,  bakery  electrochemistry,  battery  testing  and  evaluation. 

Mechanics  arid  Materials  Technology  Center:  Evaluation  and  characterization  of 
new  Ectcrials:  metals,  alloys,  ceramics,  polymers  sad  their  composites,  and  new 
forms  of  carbon;  development  and  analysis  of  thin  films  and  deposition  techniques; 
nondestructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  ccrrosica;  development  and  evaluation  of  hardened  components; 
analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  Launch 
vehicle  sad  reentry  fluid  mechanics,  heat  transfer  and  flight  dynamics;  chemical  and 
electric  propulsion;  spacecraft  structural  mechanics,  spacecraft  survivability  end 
vulnerability  assessment;  contamination,  ted  £2d  structural  control;  high 
temperature  itenaomechsaics,  gas  kinetics  and  radiation;  lubrication  and  surface 
paMaernsna. 

Spec®  end  EavfcroasBsaJ  Tcehaology  Center:  Magsetospberic,  auroral  and 
cosmic  ray  physics,  wave-particle  interactions,  magsetospheric  plasma  waves; 
atmospheric  and  ionospheric  physics,  density  and  composition  of  the  upper 
atmosphere,  remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared 
astronomy,  infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere;  effects 
of  electromagnetic  end  particulate  radiations  on  space  systems;  space  instrumentation; 
propellant  chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
atmospheric  chemical  reactions,  atmospheric  optics,  light  scattering,  state-specific 
tecsl  reactions  and  radiative  signatures  of  missile  plumes,  and  sensor  out-of-fidd- 
cf-vlew  rejection. 


